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In spite of the growing observational evidence for dark matter and dark energy in the universe, 
their physical nature is largely unknown. In fact, several authors have proposed modifications 
of Newton's law of gravity at cosmological scales to account for the apparent acceleration of the 
cosmic expansion. Inspired by such suggestions, we attempt to constrain possible deviations from 
Newtonian gravity by means of the clustering of SDSS (Sloan Digital Sky Survey) galaxies. To 
be specific, we assume a simple modification of Newton's law with an additional Yukawa-type term 
characterized by the amplitude a and the length scale A. Adopting spatially-flat universes dominated 
by cold dark matter and/or dark energy, we solve a linear perturbation equation for the growth of 
density fluctuations. In particular, we find an exact analytic solution for the Einstein - de Sitter 
case. Following the Peacock-Dodds prescription, we compute the nonlinear power spectra of mass 
fiuctuations, perform a statistical comparison with the SDSS galaxy data, and derive constraints 
in the a-\ plane; for instance, we obtain the constraints of —0.5 < a < 0.6 and —0.8 < a < 0.9 
(99.7% confidence level) for A = 5/i~^Mpc and 10/i~^Mpc, respectively. We also discuss several 
future possibilities for improving our analysis. 



PACS numbers: 04.50.-|-h 98.65.-r 98.80.Es 

I. INTRODUCTION 

Recent observations on cosmological scales suggest the 
presence of a non-zero cosmological constant or dark en- 
ergy in the universe. Adopting spatially-flat Friedmann- 
Robcrtson- Walker (FRW) models, type la supernovac [1] 
and cosmic microwave background (CMB) data [2, 3] are 
simultaneously well fit by JIa ~ 0.7 and Jim ~ 0.3, where 
f^A and r^ni denote the dimensionless cosmological con- 
stant and the density parameter of non-relativistic mat- 
ter, respectively. In fact, nearly all the available observa- 
tional data are in good agreement with the above model 
parameters. 

On the other hand, the physical origin and nature of 
such dark energy components remain to be understood. 
In this paper we will consider a possible deviation from 
Newton's law of gravity that are motivated by recent pro- 
posals for an alternative to dark energy [4, 5] or as a pos- 
sible consequence of dark energy [6-8]. While this idea 
seems exotic, the validity of Newton's law has not been 
demonstrated rigorously on cosmological scales [9, 10]. 
In addition, the presence of dark energy in the standard 
FRW model is equivalent to introducing a repulsive force 
in the cosmic expansion. In a sense, this is equivalent to 
changing Newton's law of gravity on cosmological scales 
by modifying the matter content of the universe. The 
idea behind this paper is to explore observable conse- 
quences of a modification of gravity itself. We note that 
in retrospect this is closer to Einstein's original idea in 
introducing the cosmological constant; he introduced the 
A-term on the left-hand side of the field equation. Nowa- 
days this term is usually interpreted as an additional term 



in the energy-momentum tensor on the right-hand side. 

While our current approach is entirely empirical, there 
have been several specific attempts to construct self- 
consistent cosmological models including deviations from 
Newton's law on cosmological scales. For instance, Dvali, 
Gabadadze and Porrati proposed a scenario (DGP model, 
hereafter) [4, 5] to explain the accelerating universe as 
a result of gravity leaking into extra dimensions in the 
context of a braneworld model. According to the DGP 
model, the accelerating universe is naturally explained 
without dark energy component, but rather by the modi- 
fication of Newton's law of gravity on cosmological scales. 
A variety of astrophysical and cosmological consequences 
of the DGP model were discussed by several authors [11- 
14]. Other proposals which suggest deviations from New- 
ton gravity on very large scales include a "ghost conden- 
sation" model [15, 16]. 

We would like to derive constraints on deviations from 
Newton's law that are independent of any specific models 
to the extent possible. Thus, we adopt an empirical pa- 
rameterization of the deviation from Newton's law over 
10 ^ 100 Mpc scales while the cosmic expansion is as- 
sumed to exactly follow that of the standard FRW model. 

Note that, strictly speaking, our approach is not fully 
self-consistent in that our modified gravity law has an 
effective gravitational constant Goff = Gn(1 + a) asymp- 
totically in the large-scale limit, but we still assume that 
the cosmic expansion follows the standard FRW cosmol- 
ogy with Gn- As discussed below, however, this does not 
significantly change our conclusions. 

Scalfon et al. [6] recently carried out a similar analy- 
sis in the same spirit. They derive linear theory predic- 
tions for the matter power spectrum in modified grav- 
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ity models and obtain constraints on deviations from 
the inverse-square law using 2dF (Two-degree Field) and 
SDSS (Sloan Digital Sky Survey) galaxy data. Our pa- 
per differs from theirs in the following three respects: (a) 
We do not assume a priori that any fractional deviation 
from Newton's law is small. The analysis of Sealfon et al. 
is applicable only in a pcrturbative regime (i.e., a <C 1). 
(b) Sealfon ct al. solved the linear perturbation equation 
for the growth of density fluctuations assuming an incor- 
rect scaling between wave-number and the scale factor. 
We, however, solve the original perturbation equation 
without such a scaling assumption. Finally (c) we take 
into account gravitational non-linearity by applying the 
Peacock-Dodds prescription [17]. As a result, we beheve 
that the constraints in our paper are more reliable. 

The rest of the paper is organized as follows. Section II 
describes our models for modified gravity. In Sec. Ill we 
solve the linear perturbation equation for density fluc- 
tuations in modified gravity models, and the nonlinear 
power spectra are computed in Sec. IV by applying the 
Peacock-Dodds prescription. Section V compares these 
theoretical predictions with the observed power spectrum 
of SDSS galaxies. Finally Sec. VI is devoted to conclu- 
sions and discussion. 



II. MODEL ASSUMPTIONS 

The basic equation that we solve in this paper is 
(5 + 2i7(t) j - A^<I> = 0, 



TABLE I: Parameters of cosmological models 



where 
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is a fractional mass density fluctuation. p(t) is the mean 
mass density, H{t) is the Hubble parameter, $ is the 
modified Newtonian potential and r is a proper coordi- 
nate. The dot denotes the derivative with respect to the 
cosmic time t. To solve equation (1) one has to spec- 
ify the functional form of H{t) and $. In addition, we 
will confront theoretical predictions against the observed 
power spectrum of SDSS galaxies. Thus we need an addi- 
tional assumption of the spatial biasing of SDSS galaxies 
with respect to the underlying dark matter. We adopt 
the following four major assumptions in performing the 
analysis. 

(i) For the gravitational potential we consider two 
models. The first (Model I) is to simply change the am- 
plitude of Newton's constant, Gn: 



$(r) = -GN(l + ao) / dr'^/^ 
J \r~r'\ 



(3) 



where ag is a scale-independent parameter characterizing 
the deviation. The effective gravitational constant in this 
model is simply Goff = Gn(1 + Q;o). While this simple 
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"h is the Hubble constant in units of 100km s ^ Mpc 



model is useful in understanding a basic outcome of mod- 
ified gravitational theories, it may be inconsistent with 
precise tests of Newton's gravity on small scales. There- 
fore we also consider another model (Model II) where the 
deviation is restricted to scales larger than A: 



$(r) -Gn / dr 



,3 P{r') 



1 + a 1 - e 



(4) 



Note that A is defined in proper, rather than comoving, 
length, and we consider l/i~-'^Mpc < A < 30/i~^Mpc be- 
low. Equation (4) recovers the conventional Newton po- 
tential for r <C A, and asymptotically approaches Model 
I with Goff = Gn(1 + a) for r 2> A. The transition be- 
tween these two regimes is described by the Yukawa-like 
term but our result below is insensitive to that particular 
choice. 

(ii) We assume that general relativity is valid on the 
horizon scale, i.e., that the cosmic expansion is described 
by the standard FRW model. Strictly speaking, this as- 
sumption may not be fully consistent with our modified 
gravity law on local scales (Models I and II). However, 
our primary interest in this paper is the matter clustering 
on scales of l/i~^Mpc < A < 30/i~^Mpc, which should 
not be affected by evolution of structure on much larger 
scales that are in the linear regime. In other words, even 
if we introduced another cutoff in Model II so as to re- 
cover the Newton gravity at scales ^ 100/i~^Mpc. our 
results would hardly be affected. This explains why we 
use the FRW model in describing the cosmic expansion 
while non-Newtonian gravity controls the dark matter 
clustering. For the same reason, we assume that our mod- 
ification of gravity docs not change the angular power 
spectrum of the CMB. 

In reality, we do not necessarily have to assume the 
validity of general relativity on horizon scales; the unper- 
turbed model could differ from general relativity and still 
be consistent with observations. We will defer exploring 
such possibilities for now and adopt this assumption just 
for simplicity throughout the present analysis. 

(iii) In addition, we assume that the universe is domi- 
nated by cold dark matter (CDM) and/or the cosmolog- 
ical constant. To be specific, we consider two spatially 
flat models with and without the cosmological constant: 
r^in -I- f^A = 1 (Table I). Hence we use the conventional 
models for the Hubble parameter H in equation (1) and 
the matter transfer function [18]. 

(iv) Finally, we assume that SDSS galaxies arc fair 
tracers of the total mass distribution, but allow for a 
scale-independent linear bias, 6*, for galaxies. We use 
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FIG. 1: Linear growth rate in Model I where Gcff = (1 + 
Q!o)Gn is independent of scale; (a) Sim = 1 and SIa = 0, (b) 
fim = 0.3 and = 0.7. 



assume that 6» is independent of k for simplieity. 



III. LINEAR PERTURBATION THEORY 

In this section, we focus on linear evohition of density 
fluctuation and solve Eq. (1) for the two models. Let us 
begin with Model I. In this case, the Fourier transform 
of Eq. (3) in the comoving coordinate yields 



{l\^^{x)\^ = 47rGa2p (1 + ao) 4, 



(6) 



where a is the scale factor (normalized to unity at 2 = 0) 
and (5fe denotes the Fourier transform of 5. Note that x is 
the comoving coordinate and we define k as the comoving 
wave- number. 

Substituting Eq. (6) into Eq. (1), we obtain 



2(l + "o) ( 1 



H^Sk = (7) 



in spatially flat models, where is the Hubble constant 
at the present. 

In particular, we find that Eq. (7) has an analytic 
solution in the Einstein-de Sitter model: 



_ 1 I VT+24(TT^ 
5fc cx a 4 ■» 



(8) 



Figure 1 plots the linear growth rate in Model I as a 
function of a which is independent of k. The amplitude 
in this figure is normalized 10~^ at a = 10~^. 

Next consider Model II. In Fourier space, the potential 
is written as 



the power spectrum of SDSS galaxies, Pgg(k; L^), which 
is corrected for luminosity-dependent bias and refers to 
the clustering of L* galaxies [19]. Then we attempt to 
constrain the parameter space of a and A (ao in Model 
I) by fitting to 



Pgg(fc;L,) =&*^mass(fc;a,A). 



(5) 



The scale independence of the biasing parameter is a 
fairly conventional assumption in the linear regime where 
the growth rate is scale independent. Note, however, that 
our Model II evolves in a weakly scale-dependent fashion 
even in linear theory (see Fig. 2 below). Nevertheless we 



Sk, (9) 



The linear perturbation equation for density fluctuation 
can be written as 



2HSk 
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We also find an exact analytic solution for Eq. (10) in 
the Einstein-de Sitter model: 



(5fc(a;a, A) = C*!^ 2F1 ' ^ 
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(11) 



where 2-F1 is the hypergeometrie function, and the two constants Ci, C2 are determined by initial conditions. 
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FIG. 2: Linear growth rate in Model II; (a) a-dependence for A = and k — 0.01h,Mpc~^ {Q,m = 1 and Q,a = 0), (b) 

A-dependence for a = 1 and k = O.Ol/iMpc"^ [Qm = 1 and Qa = 0), (c) fc-dependence for a = 1 and A = 10/i~^Mpc (r^m = 1 
and JIa = 0), (d) a-dependence for A = 10/i~^Mpc and k — 0.01/iMpc~^ (fim ~ 0.3 and r2A = 0.7), (e) A-dependence for a = 1 
and k = 0.01/iMpc"^ (fim = 0.3 and Ha = 0.7), (f) fc-dependence for a = 1 and A = 10/i^^Mpc (ri^ = 0.3 and Qa = 0.7). In 
all the panels, the results for a — (Newton's law of gravity) are plotted in solid curves for reference. 



The first and second terms of the right-hand side in 
Eq.(ll) correspond to the growing and decaying modes, 
respectively. 

Note that Sealfon et al. [6] introduced s = a/kX and 
attempted to solve Eq. (10) for d{k,t): 



6{k,t) ^ 5A{k,t) l + ad{k,t) 



(12) 



where (5^(fc,t) denotes the growing mode of density fluc- 



tuations for the a ~ Q case. In doing so they neglect 
higher-order terms of a and assume that c?(fc, t) can be 
written entirely as a function of s. Their scaling assump- 
tion is correct only in the Einstein-de Sitter model, but 
does not hold in general since the Hubble parameter in 
Eq. (10) depends on a and thus d cannot be expressed 
only in terms of s. Therefore we chose to solve Eq. (10) 
directly, without neglecting higher-order terms of a. 
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Figure 2 shows the Unear growth rate in Model II as 
a function of a, which is similar to Fig. 1; the left and 
right panels correspond to ilm = 1 and fljy = and to 
ilni = 0.3 and = 0.7, respectively. In this case, the 
result depends on k and A as well, and we show the a- 
, A-, and fc-dependence in the top, middle, and bottom 
panels. 

We emphasize here that A is defined in proper coor- 
dinates. Hence the gravity on the comoving scale cor- 
responding to k starts to deviate from the conventional 
Newton potential after a > kX. This is why all of the 
curves in Fig. 2 are degenerate for a < 0.01, and also ex- 
plains the fact that the linear growth rate for smaller k 
deviates earlier from the a = case in the bottom panels. 

IV. CORRECTION FOR GRAVITATIONAL 
NON-LINEAR GROWTH 

Our next task is to apply a correction for gravitational 
non-linearity to the linear predictions described in the 
previous section. The linear power spectrum in the CDM 
models extrapolated to z = (a = 1 in our definition) is 
computed according to 

PL{k;a,\) = Ar2(fc)fc"[4(a= l;a,A)]2 (13) 

where T{k) is the matter transfer function and k = |fe|. 
In the above we assume that the primordial spectrum 
obeys a power-law and fix the value of n to be unity. 
In practice we use the fitting formula [18] for T{k) 
which takes into account the baryon effect (we adopt 
51b = 0.04). It should be emphasized here that we choose 
the amplitude A so that the rms value of the top-hat 
mass fluctuation at 8/i~^Mpc equals those listed in Ta- 
ble I when a = and n ~ 1. Strictly speaking, however, 
our modified {a ^ and A ^ 0) gravity models distort 
the linear spectrum due to the factor [Sk{a = 1; a, A)]^ in 
Eq. (13). Hence the actual value of as in each model is 
slightly different from that in the a = case. This will 
be discussed in detail in Sec. V (see Fig. 6 below). 

Following the Peacock-Dodds prescription, we com- 
pute the nonlinear counterpart of Eq. (13). For that 
purpose, it is convenient to define 

A^(fc;a,A) = 47rfc3pL(fc;a,A). (14) 

Then using their nonlinear mapping [17], we transform 
(fcL, A^(fcL;a, A)) to (/cnl, A|l(^nl; a, A)), where sub- 
scripts L and NL refer to variables in linear theory and 
to those corrected for gravitational non-linear growth, re- 
spectively. In what follows, we designate and A^^ as 
the linear and non-linear power spectra. 

Consider Model II first. Figure 3 plots the linear and 
non-linear power spectra for A = 5/i~^Mpc in the A-CDM 
model with a = 0.5 and —0.5. Dashed and solid curves 
indicate Aj^{k]a,X) and A'^^{k;a, X), and A'^j^{k;a = 
0, A) is plotted in dotted curves for reference. Since all 
our models have cg ~ 1 (Table I), the non-linearity shows 
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FIG. 3: Nonlinear power spectra in Model II; (a) a — 0.5, 
(b) a = -0.5. We adopt A-CDM and A = 5ft"^Mpc for 
definiteness. Dashed and solid curves indicate the linear and 
nonlinear predictions as discussed in the text. For reference, 
nonlinear Newtonian predictions (a = 0) are also shown in 
dotted curves. 



up around wave- numbers k ^ 27r/(8/i ^Mpc). As ex- 
pected, the effect of modified gravity becomes significant 
for k < 27r/A. 

In contrast. Model I does not change the nonlinear 
power spectrum at z = relative to the ao = case if 
one adopts the identical value of erg at 2 = 0. Since the 
linear growth rate for Model I (Eq. (8)) is independent of 
k, the above statement is obvious for linear power spec- 
tra. The possible difference would come from the nonlin- 
ear evolution. If one evolves ao = and ao 7^ model 
spectra from an earlier epoch with the identical ampli- 
tude (i.e., tTg at that epoch), they lead to different values 
of (jg at z = because the ao-depcndent overall growth 
rate, and thus to different shape of the power spectra due 
to the nonlinear effect. Instead, if one chooses the iden- 
tical value of (Tg at z = 0, the power spectra of ao = 
and ao 7^ models are also identical even taking into ac- 
count the nonlinear effect. To disentangle Model I from 
the ao = model, therefore, one has to compare A^^ at 
different redshifts. Figure 4 displays the value of erg as a 
function of z. Since this value is computed by extrapo- 
lating linear theory, it can be directly obtained from the 
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FIG. 4: Linearly extrapolated values of crg(z) for different ao 
in Model I; (a) Sl^ = 1.0 and Ha = 0.0, (b) fi^ = 0.3 and 
VLk = 0.7. 

linear growth rate plotted in Fig. 1. Thus in Model I the 
structure evolves at a different rate even if the underlying 
cosmological parameters (rim and J^a) are identical. For 
example, we would need determinations of cs at multiple 
redshifts to resolve the degeneracy. 

V. COMPARISON WITH THE POWER 
SPECTRUM OF SDSS GALAXIES 

We are now in a position to explore observational 
constraints on deviations from Newton's law of grav- 
ity. In order to compare the model predictions de- 
scribed in the previous section against observations, we 
use the power spectrum of galaxies derived from 
SDSS by Tegmark et al. [19] (their Table 3 and 
http://www.hep.upcnn.edu/~max/sdss.html); they com- 
pute the power spectrum in the range 0.02/iMpc"^ < fc < 
0.3ft.Mpc~^ from a sample of 205,443 galaxies with mean 
redshift z ^ 0.1. Their measurements are in good agree- 
ment with our A-CDM model listed in Table I as long as 
a = 0. 

Figure 5 compares non-linear power spectra for Model 
II against those for the SDSS galaxies in the A-CDM 
model; Figure 5(a) shows the a-dcpendence for the A = 



FIG. 5: Comparison of power spectra between Model II pre- 
dictions and the SDSS galaxies (in the A-CDM model); (a) 
Of = 1 (dotted) and —1 (long-dashed) for A = 10/i~^Mpc, and 
(b) A = 5/i"^Mpc (long-dashed), 10/i"^Mpc (dotted), and 
30/i^^Mpc (short-dashed) for a = 1. For reference, the New- 
tonian case (q = 0) is plotted in solid curves. The data points 
with errors bars represent the results of the SDSS galaxies. 

lO/i^^Mpc case, while Fig. 5(b) illustrates that of vary- 
ing A for a = 1.0. 

In order to proceed further, we constrain a and A by 
applying Ax^ statistics [20]. Alodel II has two indepen- 
dent parameters (a and A) , and the overall normalization 
between the predictions and the data is adjusted by the 
biasing parameter 6, for galaxies which we adopt as 
a third fitting parameter. We then compute the rela- 
tive confidence levels with respect to the best-fit values 
assuming that 

Ax^(a, A) = X^(Q!, A,6,Jocal min) 

X ('^mini Aniim ^*,min) (1^) 

follows the x^distribution for 2 degrees of freedom. In 
Eq.(15), amin, Amin and 6»,min dcnotc their best-fit values 
which globally minimize the value of x^, while 6*,iocai min 
is the value that minimizes the x^ for a given set of values 
for a. and A. The results are summarized in Fig. 6; Figure 
6(a) plots the contours of Ax^(a, A); sohd, long-dashed 
and dotted lines indicate the limits at Ict, 2(t, and 3ct 
confidence levels. The amplitude of the power spectrum 
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predicted for Model II is chosen so as to lead to as = 0.9 
at z = when a 0. Thus using the linear growth rate 
plotted in Fig. 2, one can compute the contours of erg for 
Model II in the a- A plane which are shown in Fig. 6(b). 
Figure 6(c) indicates the corresponding 6*jocai min(Q:, A) 
for the SDSS galaxies with respect to each model. If the 
linear growth rate were exactly scale-independent and 
the fit were performed only over the fully linear regime, 
o'8'>*,iocai min(Q!, A) plotted in Figs.6(b) and (c) would be 
constant. In reality, however, the linear growth rate 
is slightly scale-dependent and the range which we fix, 
0.02/iMpc"^ < A; < 0.3/iMpc"\ spans a (mildly) nonhn- 
ear regime as well. Thus crsb*,ioced min(ck, A) is not exactly 
the same over the entire a-A plane. 




>.[h"''Mpc] 



VI. CONCLUSION AND DISCUSSION 

We have used a simple empirical model to parameterize 
the departure from the Newtonian law of gravity on cos- 
mological scales. We have computed the nonlinear power 
spectra of density fluctuations predicted in the model as- 
suming a spatially-flat cosmological model and derived 
the constraints on the amplitude a and the length scale 
A which characterize the additional Yukawa-type term. 
Our results are an improvement over the previous at- 
tempt by Sealfon et al. as discussed in Sees. I and III. 
For instance, we obtain the 99.7% confidence limits of 
-0.5 < a < 0.6 and -0.8 < a < 0.9 for A = 5/i"^Mpc 
and 10/i~^Mpc, respectively. 

Our Model II, Eq. (4), is intended to represent a rea- 
sonable approximation to a wide class of possible models 
which recover Newton's gravity on small scales but de- 
part from it on large scales. In contrast, Frieman and 
Gradwohl [21] and Gradwohl and Frieman [22] consid- 
ered a different non-Newtonian gravity model which devi- 
ates from Newton's law on smaller scales, but approaches 
it asymptotically on larger scales. While their specific 
model of gravity law is different from that studied here, 
they extensively discussed the effects on linear perturba- 
tion theory and explored the resulting implications for 
large-scale structure. More recently Nusser, Gubser and 
Peebles [7] considered cosmological constraints on a sim- 
ilar model of gravity which deviates on small scales. In 
order to avoid conflict with laboratory tests of Newton's 
law, they have to assume that non-Newtonian gravity 
applies only to dark matter; their model is based on the 
proposal of Ref. [8] in which the dark sector consists of 
two mutually coupled fields (dark matter and dark en- 
ergy) , and the length scale A is constant in comoving co- 
ordinate. This is why their model substantially changes 
the nonlinear clustering behavior while our model mostly 
modifies the weakly nonlinear regime. 

In either case, more specific models of non-Newtonian 
gravity are needed in order to tighten the constraints. 
Possible proposals include the DGP model in the context 
of the braneworld scenario [4, 5], and a relativistic version 
of the modified Newtonian dynamics [23]. These exam- 
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pies are attractive in the sense that they offer specific 
modified gravity laws on local scales and modified Fried- 
mann equations simultaneously in a consistent fashion. 
Therefore we could derive more stringent constraints. 

Using the SDSS galaxy power spectrum is appropriate 
for constraining the deviation around 10 ^ 100ft.~^Mpc 
scales. If a deviation from Newton's law occurs at scales 
below l/i~^Mpc, strong gravitational nonlinearity will 
be important, and nonlinear spherical collapse analysis 
and/or direct N-body simulations would be needed to 
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explore the dynamical consequences for dark halos at 
galaxy and cluster scales. If the deviation scale is much 
larger than 100/i~^Mpc, on the other hand, it is unlikely 
to leave any detectable signature on cosmic structures. 
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